Spin-polarized transport through domain wall in magnetized graphene 
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Atomically thin two-dimensional layer of honeycomb crystalline carbon known as graphene is a 
promising system for electronics. It has a point-like Fermi surface, which is very sensitive to external 
potentials. In particular, Zeeman magnetic field parallel to the graphene layer splits electron bands 
and creates fully spin-polarized and geometrically congruent circular Fermi surfaces of particle and 
hole type. In the presence of electric field, particles and holes with opposite spins drift in opposite 
direction. These phenomena are likely to be of interest for developing graphene-based spintronic 
devices. A domain wall (DW) separating regions with opposite spin polarizations is a basic element 
of such a device. Here we consider a ballistic passage of spin-polarized charge carriers through DW 
in graphene. We also discuss the analogy between the generation of spin currents in graphene and in 
relativistic quark-gluon plasma, where the spin-polarized current is responsible for the phenomenon 
of charge separation studied recently at RHIC. 

PACS numbers: 73.63.-nb, 73.40.-c, 72.25.-b, 75.25.Mk 



I. INTRODUCTION 

The remarkable properties of graphene^ at present 
attract attention of many researchers. Its honeycomb 
two-dimensional (2D) crystalline order is extremely ro- 
bust. In view of the well-known Peierls-Landau argument 
proving the thermodynamical instability of isolated 2D 
crystals^, robustness of graphene monolayers may look 
somewhat surprising. While graphene sheets in natu- 
ral and synthetic graphite materials are not isolated but 
are supported on 3D substrates^^^, binding of the mono- 
layers is so weak that they can be easily exfoliated and 
in many cases appear approximately isolated from the 
substrate^i^i^il. 

Adding to the enthusiasm were recent discoveries 
of exceptionally high electronic quality of graphene 
monolayers''^. Carefully prepared samples show am- 
bipolar electric field effect with carrier mobilities ex- 
ceeding 10^ cm^/V/s for electron/hole concentrations up 
to ~ 10^^ cm~^ .— High quality graphene samples at 
present can be obtained either by using graphite exfolia- 
tion, which results in graphene pieces with 1 to 100 jiia 
linear dimension^, or by epitaxial growth on SiC via sil- 
icon sublimation, which yields macroscopic mosaic layer 
with micron-size crystalline domains^^. 

High electron mobility implies ballistic charge trans- 
port and electronic phase coherence on the micron length 
scale, which is comparable with high-quality semiconduc- 
tor heterostructures traditionally used for the studies of 
the quantum Hall effect (QHE)^^. Moreover, charge mo- 
bilities in graphene are only weakly temperature depen- 
dent, being probably limited by sample imperfections and 
size effects even at room temperature^. Hence, not only 
QHE was indeed observed in graphen o^^i^^ , but it was 
also found to survive up to 300 K^^, indicating that the 
electrons in graphene form a quantum gas even at room 
temperature. 

The exceptional electronic properties of graphene such 



as high charge carrier mobility, long mean free path 
and coherence length and ability to support high cur- 
rent densities, exceeding lO' A/cm^,' make it a promis- 
ing candidate for nano-scale electronics. In addition, 
its unique electronic structure provides a playground for 
the studies of (2+l)-dimensional (space+time) quantum 
electrodynamics. The crossing of the energy bands as- 
sociated with two different sub-lattices A and B of the 
graphene's honeycomb crystal lattice results in the en- 
ergy spectrum of electron and hole quasi-particles, which 
is linear in momentum hk [see FiglTJa)], e{k) = vphk, 
where vf — 10^ m/s is the Fermi velocity. This has been 
observed experimentally in transport measurement □ ^ ^ i "'^ ^ 
and in angle-resolved photoemission-. 





FIG. 1: (a) Electronic band structure resulting from the sp^ 
C-C bonding in the hexagonal carbon layer of graphene. In 
zero magnetic field the filled tt and the empty tt* bands 
meet at a single point, resulting in a linear 2D dispersion, 
s(k.) = vphk, characteristic of 2D relativistic Dirac fermions 
(vf is Fermi velocity), a and b define triangular lattice of 
a honeycomb graphene layer containing two C atoms. Fermi 
"surface" consists of two inequivalent points K and K' (val- 
leys), (b) Magnetic field H parallel to graphene layer intro- 
duces Zeeman splitting q/jbH between the bands with parallel 
(P) and antiparallel (AP) spin. P and AP bands acquire con- 
gruent Fermi-surfaces of hole- and electron-type respectively, 
whose radius is hkF = qiibH / {2vf)- 



A gapless linear 2D spectrum of electron and hole 
quasi-particles belonging to two sublattices implies that 
charge carriers in graphene can be formally described as 
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two-dimensional relativistic chiral fermions with spin and 
with pseudo-spin accounting for the two-sublattice band 
structure. A straightforward consequence is the conser- 
vation of the pseudo-spin chirality of the quasiparticles 
defined as the pseudo-spin projection on the momen- 
tum, p • T. The orbital motion in magnetotransport and 
QHE experiments can be described by the "truncated" 
2D Dirac equation for massless fermionsi^!^!^, 

(p + -A)Tlljr' = eipr' , (1) 

c 

where = (r^^ry)^ are Pauli matrices acting in the 
pseudo-spin space. They account for the two-sublattice 
nature of graphene's honeycomb lattice and the resulting 
composite structure of the dispersion cone around each of 
the Fermi points, ipr' is a rank two spinor wavefunction 
and p is the 2D momentum operator^^^i^iii. 

Eq. ([1]) assumes degeneracy with respect to spin and 
valley indices, which are taken into account simply by 
multiplying the number of states by 4. Spin degener- 
acy is usually justified by the fact that typical Zeeman 
electronic level splitting induced by the laboratory mag- 
netic field is indeed very small. However, when magnetic 
field is applied parallel to the plane of graphene, orbital 
motion and Landau quantization are irrelevant and it is 
the lifting of spin and valley degeneracies that becomes 
important^^. Such situation is in fact of immediate in- 
terest for possible spintronic applications, which employ 
the spin degree of freedom of Dirac fermions in graphene. 
With recent measurements showing spin coherence scale 
in graphene exceeding 1/im^^, there is a growing recogni- 
tion of the potential of this appro achiSiSSiSi. Moreover, 
it can be envisioned that spin-dependent splitting of the 
electronic levels in graphene could be induced by an effec- 
tive magnetic field resulting from magnetic proximity ef- 
fect in graphene in contact with a ferro/antiferromagnetic 
substrate^^'^^. Such "exchange" field acts only in the spin 
sector and can be much stronger than magnetic fields 
available in the laboratory, inducing level splitting suffi- 
cient for room-temperature device applications. 

Here we consider transport of chiral Dirac fermions in 
graphene in the presence of such an in-plane magnetic 
field and their passage through a boundary between two 
regions with different field orientations. The latter might 
be associated with a domain wall in the magnetic layer 
of graphene-magnet (GM) heterostructure, and presents 
a basic element for spintronic applications. 

II. EFFECT OF PARALLEL MAGNETIC FIELD 

When external magnetic field is applied to graphene, 
its parallel component acts only on spin degree of free- 
dom, while the perpendicular component couples both to 
spin and to the orbital motion as described by Eq. ([T]). 
Hence, the action of the parallel field is equivalent to 
band splitting by the "exchange field" arising from mag- 
netic proximity effect induced by the magnetic substrate. 



Such proximity-induced field, regardless of its direction, 
does not couple to the orbital motion. 

Point-like Fermi surface makes electronic properties of 
graphene extremely sensitive to external potentials. Ex- 
periments show that application of a moderate gate volt- 
age results in an appearance of finite charge carrier den- 
sity, which is proportional to the magnitude of applied 
electric field. The type of these induced carriers, revealed 
by the sign of the Hall effect, depends on the polarity of 
the gate voltage^' This can be easily visualized by 
considering the electronic spectrum in graphene shown 
in Fig. 1(a), where the Fermi level is shifted by applying 
the gate voltage, thus inducing a circular Fermi-surface 
of particle or hole type, depending on the polarity. 

External magnetic field H splits electronic band struc- 
ture in graphene according to spin. Chemical potential 
for one spin polarization is increased by the amount equal 
to Zeeman energy ^/i^i^/2, while for the other it is de- 
creased by the same amount. Fig. [ifb) 2 is the 
spectroscopic Lande factor for electrons in graphene, (Ib 
is Bohr magneton). As a result, there appear identical 
circular Fermi surfaces, of particle type for spin antipar- 
allel to magnetic field, and of hole type for spin parallel 
to it. The radius of these Fermi surfaces, hkn^ is propor- 
tional to the magnetic field, 

hkH=g^BH/{2vF). (2) 

The difference in filling of the two spin states results in 
small Pauli paramagnetic moment and total charge car- 
rier density at the Fermi level, 

n{H) = {g^BHf/{2nh\l). (3) 

In order to achieve non- negligible carrier densities, 
magnetic fields yielding Zeeman splitting of hundreds 
of Kelvins or more are required. While such fields 
can not be produced by solenoids, they might be in- 
duced by magnetic proximity effect in graphene-magnet 
multilayers2^i2^. 

III. THE EFFECTIVE HAMILTONIAN 

In this section we review the low-energy Hamiltonian 
describing electrons in graphene in the presence of Zee- 
man field. Recall that the Fermi surface of undoped 
graphene contains two non-eqiuvalent points, K and K' ^ 
giving rise to a valley degeneracy. At each of these 
points the wave-function is a pseudo-spinor in the two- 
dimensional space of A and B sublattices and a spinor in 
the spin angular momentum space. We use Pauli matri- 
ces Tx^y^zi and cFx^y^z to refer to the sublattice pseudo-spin 
and the "usual" spin, respectively. With these notations, 
the effective Hamiltonian takes the form 

^ = vfP-t + B(x)o-, (4) 

where B = ^/iBH(x)/2, and spatially varying magnetic 
field H(x) couples to spin degree of freedom. As it was 
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discussed above, we envision that this magnetic field can 
be induced by the proximity effect in graphene due to 
the super-exchange interaction with the magnetic layer 
contacting the graphene sheet^^^^^. As a consequence, 
magnetic field considered in Eq. (|4]), irrespective of its 
direction, acts only on spin degree of freedom of quasipar- 
ticles in graphene. For the case of spatially homogeneous 
magnetization of graphene-magnet heterostructure, the 
spin component directed along the effective field is con- 
served. Hence, spin-polarized carriers maintain their po- 
larization. Such structure can be used to transport spin- 
polarized currents. The basic element allowing manipu- 
lation of spin-polarized currents in graphene-magnet het- 
erostructures is the region where magnetic field changes 
its direction, namely the domain wall (DW). In the next 
section we analyze transmission of the spin-polarized car- 
riers through the domain wall. 



a) 




FIG. 2: Two regimes of spin-polarized electron transport 
through a domain wall between the regions with opposite spin 
polarizations in magnetized graphene. (a) Passage through 
narrow (< 10 nm) domain wall is mediated by Klein tun- 
neling of electron-hole pairs, (b) Transport through thick 
(> 100 nm) domain wall allows electron's magnetic moment 
to follow the magnetic field direction. 



IV. PASSAGE OF SPIN-POLARIZED DIRAC 
FERMIONS THROUGH A DOMAIN WALL 

In order to understand transport properties of an in- 
homogeneously magnetized graphene heterostructure, we 
analyze ballistic passage of spin-polarized carrier through 
the lateral domain wall separating two regions of opposite 
magnetization. Fig. [2l The transmission through the do- 



main wall is characterized by the amplitudes of spin flip 
and non-spin-flip processes. These amplitudes determine 
spin-polarized transport through DW and their knowl- 
edge is important for devising spintronics applications of 
the considered heterostructure. 

For definiteness, let us consider electrons in the pres- 
ence of the magnetic field B pointing along and opposite 
to z-axis on different sides of the domain wall located in 
the stripe region xq < x < xq -\- L. We further assume it 
to rotate uniformly within the DW (see Fig. [2|). Specif- 
ically, we represent the magnetic field B = Bn{x) with 
the unit vector 

{-Uz X < Xo 

—Uz COS + Uy sin 'd{x) xq < x < xq -\- L 
Uz X > Xq ^ L 

(5) 

describing its rotation within the DW. In Eq. ([5]), 
stands for the unit vector pointing in the i-th direction, 
and the rotation angle is taken to be linear in the lateral 
coordinate x, = 7r{x — xo)/L. 

We begin with the qualitative discussion of the car- 
rier passage in the case of normal incidence. Due to the 
conservation of (pseudo-spin) chirality in the Klein tun- 
neling phenomenoni^ the backscattering is absent in this 
case. The carrier passes the DW during the passage time 
tp = L/vf^ where L is the characteristic width of the 
DW. Within the DW region the spin of the carrier ex- 
periences the time dependent torque and undergoes the 
Larmor precession. The parameter controlling the spin 
dynamics is the ratio of the passage time to the spin pre- 
cession period, 

r] = 2BL/7TVF. (6) 

For thin DW, 77 <C 1, and in the absence of the interval- 
ley scttering, Xds/L <C 1, where XdB is the carrier's de 
Broglie wavelength, the spin-flip probability is small, cor- 
responding to a small precession angle. In the opposite 
limit of thick DW, 77 1, the spin follows the varying 
Zeeman field inside the DW adiabatically and the non- 
spin- flip probability is small. As a result, the carrier 
preserves its alignment with the field and reverses its po- 
larization upon passing through the DW. Although the 
scattering for the arbitrary angle is complicated by the fi- 
nite backscattering amplitude, the basic physical picture 
presented above still holds and allows us to construct a 
scattering theory in the general case. 

To solve the scattering problem we construct the scat- 
tering state, 

// X i^ii^^z) ^^aa'^a'i^^^)^ X < Xq 

where '^^(x, z) with s = i^r^t denotes incoming, reflected 
and transmitted waves, and the subscript a = ± refers to 
the spin up (down) polarizations, respectively. Due to the 
translational symmetry in the z-direction, the scattering 
state in Eq. ([7j) is the eigenstate of the 2;-component of the 
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momentum and can be labeled by its eigenvalue Pz , mak- 
ing our problem effectively one-dimensional. We present 
wave functions entering Eq. ([7j) in the following form 



s = z, r, L (8) 



Here Xa is the spin wavefunction, i. e. the spinor satis- 
fying azXa = (^Xa^ and (Pa,pz ^ pseudo-spinor in the 
sublattice space. 

Capitalizing on the particle-hole symmetry of the prob- 
lem, in what follows we only consider the case of inci- 
dent quasiparticles with ^ > 0. The incoming wave is 
an eigenstate of the Hamiltonian ^ for x < xq. At a 
fixed energy, the majority (spin down) and minority (spin 
up) carriers have Fermi momenta p-^ = {E -\- B)/vf and 
P- = \AE\ /vf, respectively. Here AE = E — B can be 
both positive and negative, the latter case corresponds to 
the hole-like quasiparticles. The kinematic constraint for 
the incoming spin up (down) electrons reads \pz\ < p^. 
Introducing the notation 



u±{Pz) = V(l+P./P±)/2 



(9a) 
(9b) 



for the pseudo-spinor (pp^ = (cos 7(^^)72, sin 7(^^)72) 
forming angle ^{pz) with z-axis, we write the pseudo- 
spinor of the incoming electron as 



U-{Pz) 
V-{Pz) 



-AE 



V-{Pz) 
U-{Pz) 



(10a) 



where Oae and 0-ae are step functions distinguishing 
cases of particle- and hole-like carriers, respectively. In a 
similar fashion we write for the refiected wave 



r 



yAE 



yAE 



U-{Pz) 


+ 0-AE 


' V-{Pz) 


_ -V-{Pz) 


. -Ul{pz) _ 


U+{Pz) 


+ 0-AE 


'uliPz) 


_ -V+{pz) _ 


_-yUPz) _ 



. (10b) 



The square roots in Eq. ([9]) are defined as having posi- 
tive imaginary part for negative argument. This choice, 
together with the sign and conjugation convention in 
Eq. (|10bp . insures that for p_ < \pz\ < p+ the wave 
function of the minority (spin up) carriers decays expo- 
nentially away from the DW, namely, it is an evanescent 
wave. The transmitted waves are given by 



U+{Pz) 



^AE 



U-{Pz) 
V-{Pz) 



-AE 



vl{Pz) 
uliPz) 



(lOc) 



Equations (p!Q|) when substituted in ([8]) give the explicit 
expressions for the incoming, refiected and transmitted 
spinors in the most general scattering state of Eq. ([7j) . 



In order to find the transmission and refiection ampli- 
tudes, the transfer matrix T matching the wave function 
at the two boundaries of the DW 



^{xo)=T^{xo^L) 



(11) 



has to be found. To this end, we solve the Dirac equation 
inside the wall, 

{—ivpTxdx + B(Jz cos^{x) — B(jy sm'd{x)) ip 

= (E-VFPzTz)^. (12) 

with the initial condition specifying the wave function 
at X = xq. The equation ([12]) is formally equivalent to 
Rabi problem of spin coupled to the oscillating magnetic 
field^^ with coordinate x playing the role of the time. The 
field is turned on at the "time" x = xq and turned off 
at "time" x = xq -\- L. Exploiting this analogy we solve 
Eq. (p!2|) by the transformation to the rotating reference 
frame, 

tlj{x) = exp (^—i'd{x)^^ 'ij){x) (13) 

such that the field seen by the transformed spin is sta- 
tionary. Substitution of Eq. (p!3|) into Eq. ([12]) gives 



TTVf 

~2L 



Bgz) i^{x) 

= {E-VFPzTz)i^{x). (14) 



We notice that static magnetic field now appears effec- 
tively as an operator in the pseudo-spin space. We can 
rewrite Eq. ([T4|) in the form 



(15) 



where ( — {x — Xq) / L. The four-by-four matrix on the 
right hand side of Eq. ([T5]) reads 



-Gx + r]TxCFz - r]Tx{e - CzTz) 



(16) 



with dimensionless energies e = E/B and = vfPz/B 
and the parameter 77 defined in Eq. ([6]). The formal so- 
lution of equation ([T5|) is 

ij{xQ + L) = exp V^(^o) . (17) 

Combining equations ([TT]) . ([T3|) and ([T7|) we obtain 

f = exp(z|i)i(7^. (18) 

The exponentiation in Eq. ([T8|) can be performed explic- 
itly as follows 



cos 



lA 



sin (f Vc± A) 



Vc±A 



(19) 
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where notations 



(20) 



A = 2r?\/e2(l + 7?2) 



have been introduced, and 



P± 



2A 



(21) 



are projection operators onto the subspaces of the eigen- 
values ±A of the matrix 



M 



i2 _ 1 Tri^ 1 Tri^ 
4 4 



(22) 



Equations ([IHl), ([T9l), ([2T]), and ([22|) give the transfer ma- 
trix T introduced in equation (pT]) as a third degree poly- 
nomial in matrix A. For an arbitrary angle of incidence 
the transmission and reflection coefficients are found by 
imposing the matching condition in Eq. ([TTj) on the scat- 
tering state of Eq. ([7j). The resulting system of four linear 
equations determining four coefficients rQ,Q,/, to^o^' can be 
easily solved. For the incoming (spin up) minority car- 
riers the probabilities of the passage with and without 
spin flip are given by 



{Pz/P-f 



(23) 



Although the solution of linear matching equations is 
straightforward, the final expressions are somewhat cum- 
bersome and below we present the results graphically. 
The transmission probabilities for the incoming minority 
(spin up) and majority (spin down) species, respectively, 
are shown in Figs. [3] and [H In both figures, panels (b) 
and (d) present the transmission with the spin flip, and 
panels (a) and (c) the transmission without the spin flip. 
Panels in the upper row show the case of the thin DW, 
L = 5, 7^ <C 1, and panels in the lower row the case of the 
thick DW, L = 50,7?> 1. 

Unlike the case of the normal incidence, for the inci- 
dence at an arbitrary angle the chirality is not conserved, 
leading to a finite backscattering probability. Hence, par- 
ticles passing the DW experience spin-dependent reflec- 
tion and refraction. The Snell's law relating the angles of 
propagation of the incoming and the outgoing particles 
to the refraction indices of the two media reads 



n^{x < xq) sin^^(x > xq ^ L) 



n^{x > xo -\- L) sm6^{x < xq) 



(24) 



The conservation of the z-component of the momentum 
gives for the ratio of the refractive indexes. 



n^(x < Xo) B 
n^{x > xo^ L) ~ E±B ' 



(25) 
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FIG. 3: Transmission probabilities for the incoming minority 
(spin up) carriers for two different DW thicknesses, L. The 
probability is shown on a linear gray scale from (dark) to 
1 (white). The parameters used are B — 0.05 and vf — 1. 
Panels present the probability of transmission (a) without the 
spin flip, L = 5 ; (b) with the spin flip, L = 5; (c) without 
the spin flip, L = 50; (d) with the spin flip, L = 50. The 
parameter 77 defined by Eq. (|6]) is 77 = 0.16 for the thin DW, 
(a), (b), and ry = 1.6 for the thick DW, (c), (d). 
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FIG. 4: Transmission probabilities for the incoming majority 
(spin down) carriers for the set of parameters used in Fig. [3] 



The behavior of this ratio for two spin polarizations as a 
function of energy is shown in Fig. [5l 

It follows from Eq. ([25|) that the DW has a different 
effect on the spin minority and the spin majority carriers. 
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FIG. 5: Ratio of the refraction indices on the two sides of DW 
for the majority (-) and the minority (+) incident particles. 



The minority carriers passing the DW experience the in- 
crease of the optical density, while the majority carriers 
experience the decrease of it. An interesting regime oc- 
curs when the energy of incoming particles E < B. In 
this case the ratio in Eq. ([25]) becomes negative and re- 
sults in the spin dependent Veselago lens effect, similar 
to that discussed in Ref. 26. 

The above spin-optics arguments are useful in under- 
standing the results shown in Figs. [3] and [H For the in- 



cidence at a shallow angle the probability of non-spin-fiip 
passage is suppressed as particle trajectory even in the 
narrow DW becomes long. The difference between the in- 
coming minority and majority species in Figs.[3l(a,b) and 
|4](a,b) results from the different refraction coefficients for 
the two species imposed by the kinematical constraints. 
The refraction coefficient ratio for minority carriers is 
n+ < 1, and their trajectories bend so that the path in- 
side the DW shortens. For the majority carriers, on the 
other hand, n_ < 1, and the trajectory bending leads 
to the longer path inside the DW. Therefore, the effect 
of magnetic field inside the DW and the probability of 
transition without the spin flip are enhanced for the mi- 
nority carriers, Fig.[3l(a), and suppressed for the majority 
carriers. Fig. [H^a). 

The Fabry-Perot pattern of transmission seen in 
Figs, ^d) and W^d) is a consequence of an interference 
of multiple reflections inside the thick DW. 



A. Normal Incidence 

Our results are substantially simplified in the case of 
normal incidence, when the momentum component Pz 
vanishes. With = the equation (p!9|) reduces to 
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The last equation gives for the transfer matrix 

f = ia.cos^vW- ^+!!^^ sin^vW, (27) 

where the overall phase factor e*^^^/^ has been omitted. 
In the present section we focus on the transmission prob- 
abilities. It has to be stressed however that the phase of 
the transition amplitude is also of interest, especially if 
the magnetization vector B completes one or more rota- 
tion circles inside the DW. Under the conditions of adi- 
abatic spin transfer, this phase is geometric^^, see the 
discussion of geometric Berry phase for the DW passage 
in App. El 

In the case of normal incidence the chirality is a good 
quantum number, [T^Tx]- = 0. This ensures the ab- 
sence of backscattering (Klein tunneling phenomenon). 
For the incoming particles with £^ > we have Tx = +1. 
Therefore, dynamics in the case of normal incidence oc- 
curs in the spin sector only. The particle traveling inside 
the DW experiences the action of magnetic field rotating 
with the frequency co = ttvf/L. The probability of pass- 
ing the DW without the spin flip is given by the diagonal 



element of the transfer matrix 

T++ = (-^) sin^ {ntp/2) . (28) 



Here, Qtp = iry/l^r]'^ is the rotation angle accu- 
mulate d by the spin pre cessing at the Rabi frequency, 
Q = + during the passage time tp = 

L/vp. It follows from Eq. ([28]) that the polarization of 
the impinging particle is not influenced by the DW in the 
case of the thin wall, 77 ^ 1. In the opposite limit of the 
thick DW, 7? ^ 1, electron spin adjusts adiabatically fol- 
lowing the direction of the magnetic field slowly varying 
inside the DW. 

Our results for the case of normal incidence are in 
agreement with Ref. [lH, where neutron polarization 
change in the course of the passage through the ferro- 
magnetic domain wall is analyzed. We note that in the 
non-relativistic case^^ the absence of backscattering is 
an approximation, which is valid for neutrons with high 
enough energy. This approximation fails for low energy 
particles, as inside the DW particles experience a force 
which results from the magnetic field gradient inside the 
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DW. In the present, relativistic, case however the carri- 
ers with arbitrarily low energy can not be reflected by 
the DW because of the chirality conservation. Therefore, 
in the relativistic case the thick domain wall is 100% ef- 
ficient in flipping spins of particles incident at 90°. 

Until now we have discussed the case of the DW 
with well-defined, abrupt boundaries, where the region 
of magnetic field variation is limited to a finite interval 
[see Eq. (|5])]. While in most cases this is a reasonable 
description (in particular, for patterned structures), in 
some experimental realizations of spintronic devices the 
boundaries of the DW may be smooth and not well de- 
fined. To clarify the significance of the above distinc- 
tion in the DW structure we consider the DW with the 
magnetic field direction n = — Uz cos i9{x) + Uy smi^{x) 
with the angle i^{x) following the Rosen- Zener profile, 
dx^{x) = (tt/L)/ cosh(7rx/L). In this case the non-spin- 
flip transmission amplitude can be found exactly^^. 



sech^ 



(29) 



It follows from comparison of Eq. ([28|) and Eq. ([29|) that 
the DW with smooth boundaries polarizes the incoming 
carriers even more efficiently than the DW with abrupt 
boundaries. Hence, we argue that both for thin DW, 
77 <C 1, and thick DW, rj ^ 1^ cases our conclusions are 
valid for the DW of an arbitrary shape. 
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FIG. 6: (Color online) The partial conductances 
per unit sample width, W, as a func- 
tion of the chemical potential, /x, for DW of two different 
lengths L and spin up {(3 = +) and spin down ((3 = — ) 
incoming carriers, B = 0.05. (a) L = 5, /3 +, (b) L = 50, 
/3 = +, (c) L = 5, /3 = — and (d) L = 50, /3 = — . In all panels 
the off diagonal conductances with a ^ [3 corresponding to 
a spin flip processes are shown by dashed (blue) line, the 
diagonal conductances are shown by solid (red) line. The 
thin (black) line shows the conductance in the absence of the 
DW per one spin and one valley as a reference line. 



V. CONDUCTANCE IN THE BALLISTIC 
TRANSPORT REGIME 

In the high quality graphene devices the mean free path 
is comparable to the characteristic sample size. Under 
such conditions the transport through DW structures is 
ballistic. The conductance can be obtained within the 
Landauer approach. In spintronic devices we consider 
the spin selective transport, as they manipulate the cur- 
rents of the electrons of different polarizations indepen- 
dently. Two terminal conductance is given by the sum of 
the transmission probabilities over all active conductance 
channels. Both majority and minority spin channels give 
rise to currents of carriers of both polarizations. We in- 
troduce the spin dependent conductance Gcx(3 to denote 
the contribution of incoming carriers with spin (3 in the 
source channels to the current of carriers with spin a in 
the drain channels, 

Ga/? = j ^ ^T^piP.) , (30) 

where we have taken into account the 2-fold valley de- 
generacy. The above definition is meaningful due to the 
conservation of z-component of spin away from the DW. 

Partial conductances Ga/s/i'^e'^W/h) obtained from 
Eq. (|3Q|) are plotted in Fig. [6] as a function of the chemi- 
cal potential for two thicknesses of the DW, L = 5, (a,c) 
and L = 50, (b,d). The partial conductance for the pos- 



itive spin polarization in the source, = +, is shown in 
the upper row. Fig. [6^a,b). That for /3 = — is shown in 
the bottom row. Fig. [6l(c,d). Solid lines corresponds to 
the spin flip processes, namely a ^ 13 and dashed lines 
represent the diagonal conductances with a = (5. 

The common feature of the curves shown in Fig. [6] is 
the conductance growth when the chemical potential /i 
exceeds half of the spin splitting, B. This is clearly due 
to the increase of the number of conducting channels. 
Secondly, the increased thickness of the DW stimulates 
spin flip processes. This is a consequence of the adia- 
batic transfer of the spin inside the thick DW discussed 
in Sec. IV. In addition, the only conductance not van- 
ishing at the special point ji = B in the case of the thick 
DW is G+-, see Fig. [6fd), dashed line. This is eas- 
ily understood from the following consideration. For the 
thin DW the spin is approximately conserved. For that 
reason, the spin minority (majority) carriers have van- 
ishingly small number of incoming (outgoing) channels 
at /i ~ 5, leading to a small conductance. In the case 
of the thick DW the spin majority carrier can remain 
spin majority carrier by adiabatically adjusting (revers- 
ing) its spin polarization. This yields finite conductance 
dit fi ^ B. The specifics of the point ji = B described 
above makes hetero-structures with the Dirac spectrum 
of quasi-particles promising candidates for spin manipu- 
lation of the currents in spintronics. 
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VI. SUMMARY AND DISCUSSION 

In the present paper we have analyzed the passage 
of spin-polarized Dirac charge carriers through the DW 
in graphene-magnet hetero-structure. We have calcu- 
lated the transmission and the reflection probabilities as 
a function of the energy of the incoming particles, which 
is determined by the average chemical potential /i in the 
graphene sample, for the DW of different thickness. The 
knowledge of the transmission amplitudes has allowed us 
to calculate the conductances of different spin channels 
in two-terminal geometry. The spin-polarized transport 
depends crucially on the thickness of the DW. Below we 
discuss the main features of our results and their poten- 
tial applications in graphene-based spintronic elements. 

We have considered two limiting cases of thin and thick 
DW. In the case of thin DW the spin dynamics inside the 
DW only occurs for shallow incidence angles. Aside from 
this special case, the spin is approximately conserved and 
the transmission is governed entirely by the kinematics 
of the relativistic Dirac quasiparticles in graphene, estab- 
lishing direct connection with the problem of Klein tun- 
neling and chiral dynamics in 2D quantum electrodyman- 
ics. A similar problem has recently been considered in the 
context of p — n junctions in graphene devices^^^^^^^^. In 
our case the DW presents a p—n junction for the majority 
and Si n — p junction for the minority carriers. The non- 
spin-flip transmission for E < B is allowed through the 
particle-hole transmutation - the Klein tunneling phe- 
nomenon. 

The transmission properties of the thin DW can be 
nicely understood in the context of spin-optics^. As it 
follows from the ratios of the refraction indices in Fig. O 
the trajectories of the minority (spin up) carriers bend 
inward, while those of the majority carriers bend out- 
ward. This difference is most pronounced near E = 
where the refraction index changes in sign, becoming neg- 
ative for both spin polarizations dit E < B. Near this 
point the transmission of the majority carriers through 
narrow DW vanishes and the corresponding refraction 
index diverges. The index of refraction for the minor- 
ity carriers, on the other hand, is close to 0, as they 
can only pass through the DW near the forward direc- 
tion. Hence, there is a giant birefringence of carriers with 
different spin polarizations, which could be employed in 
spin-selective transport devices. At E ^ B, the spin ma- 
jority carriers undergo the total internal reflection due to 
the kinematic constraint, which occurs in a wide angu- 
lar interval corresponding to dark areas in Fig. [4](a). In 
the regime close to the total reflection the length of the 
particle's trajectory even inside a narrow DW becomes 
increasingly large. Hence, the probability of the spin- flip 
transmission becomes significant, see Fig. HJ^b). In this 
regime the angular aperture of the spin minority carriers 
becomes small, and thin DW is an efficient spin polarizer. 

In the case of the thick DW the passage is governed by 
the spin dynamics inside the domain wall and the con- 
ductance is controlled by the adiabatic nature of the spin 



transport. Independent of their energy, both majority 
and minority carriers simply flip their spins, remaining 
in their respective channels, so that thick DW acts as a 
spin-flipper for spin-polarized currents. As in the case 
of the thin DW considered above, the effect of the DW 
on the transport is most pronounced when the chemi- 
cal potential is tuned to near the spin splitting, ji ^ B. 
In this case the Fermi surface of the minority spin car- 
riers in the region x > xq ^ L shrinks to a point and 
the current is carried by the majority species. In con- 
trast to the case of the thin DW, the majority carriers 
in the region x < xq are transformed to the majority 
carriers in the region x > xq ^ L hy adiabatically adjust- 
ing their spin. This gives finite conductance ioi ji ^ B^ 
which is off-diagonal in spin, see Fig.[6^d). Thus, a spin- 
transistor action could be achieved by virtue of adjusting 
the chemical potential in a gated device with magnetic 
layer coupled to the graphene layer. 

Finally, manipulating the spin polarization of elec- 
trons in graphene by the means of Zeeman band splitting 
such as discussed in this paper depends crucially on the 
strength of the polarizing magnetic field. In the case of 
the passage through a DW, this strength also determines 
the relevant physical thickness distinguishing the cases of 
thick and thin DW, Eq. (|6]). In the case of the labora- 
tory magnetic field of Hb ^ IT created by a solenoid, 
the Zeeman band splitting is B ^ 0.1 meV. The condi- 
tion ?7 ^ 1 requires extremely thick DW, L* ^ 20/im. 
Moreover, such field results in a negligible spin-splitting, 
which corresponds to a temperature of only ^ 1 K, and 
negligible charge carrier density, uh ^ 10^ cm~^, Eq (j3|). 

However in the case of a spin-dependent band splitting 
induced by the magnetic proximity effect in graphene- 
magnet hetero-structure the corresponding effective spin- 
polarizing field could be expected to be of the order 
of hundreds, or even a thousand Kelvins. Physically, 
it could be estimated from the characteristic ordering 
temperature of the magnetic layer, and corresponds to 
10 meV to 100 meV spin-dependent band splitting. 
Then, the effective "exchange" field could be as large as 
~ 10^ T, resulting in uh ^ 10^^ cm~^ and L* ~ 200 nm. 

Is it reasonable to expect such large magnetic proxim- 
ity effects in graphene-magnet heterostructures? Theo- 
retical estimates for the graphene layer in contact with 
EuO, one of the few insulating Heisenberg ferromagnets, 
predict band splitting of ~ 5 meV— . This agrees well 
with the Curie temperature Tc = 69.3 K for this mate- 
rial. For room-temperature magnetic materials, the split- 
ting could be expected to be proportionally higher. For 
example, one can envision using NiO [111] film in con- 
tact with graphene. NiO is an antiferromagnet with the 
Neel temperature Tn ^ 500 K, where ferromagnetically 
aligned Ni layers alternate in the [111] direction. In the 
artificial [111] nano-layer CoO/NiO and NiO/Fe304 su- 
perlattices the proximity effect-induced increases of the 
ordering temperature by hundreds of Kelvins have been 
observed^^. It does not seem unreasonable to extrapo- 
late this effect to graphene-magnet heterostructures. Fi- 
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nally, the technology of manufacturing graphene-based 
heterostructures is developing quite rapidly. The growth 
of the atomically smooth epitaxial MgO film on graphene 
has recently been reported^^, as were nonvolatile mem- 
ory devices obtained by covering graphene with a fer- 
roelectric layer Therefore, manufacturing graphene- 
magnet heterostructures in order to achieve manipula- 
tion with spin-polarized currents such as considered in 
this paper looks reasonable and promising approach to 
be attempted experimentally. 

In closing, we would like to note an interesting analogy 
between the generation of spin-polarized electric current 
in graphene hetero-structures and the separation of elec- 
tric charge of quarks in strong external magnetic fields in 
the presence of the topological charge in the quark-gluon 
plasma produced in relativistic heavy ion collisions'^. 
There is a recent evidence'^ for this effect from STAR 
Collaboration at RHIC (BNL). 
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APPENDIX A: GEOMETRICAL PHASE OF 
CARRIERS PASSING THE DW 

The purpose of this appendix is to illustrate the gen- 
eral concept of geometrical phases^^ for the spin carriers 
passing the DW described by the Hamiltonian of Eq. (|4]). 
We consider the normal incidence case where the dynam- 
ics occurs in spin sector only. Here we are interested in 
the limit of adiabatic spin transfer realized in sufficiently 
thick DW. The direction of the magnetic field is a pa- 
rameter of the system changing slowly inside the DW, 



B = B(x) + 5^n^, 
where B(x) = Bn{x), and 



(Al) 



Iz X < Xo 

n{x) = ^ Uz cosi3{x) — Uy sini3{x) xq < x < xq -\- L . 



X > Xo 



(A2) 

In the last equation we assume i9{x) = AO{x — xo)/L. 
For = 27r, the magnetic field defined by Eqs. ()Aip 
and ()A2p acting on a particle passing the DW completes 
a circle, sweeping the conical surface in B-space. It is 
identical for x < xq and x > xq -\- L. Therefore, the 
wave function adiabatically following the instantaneous 
eigenstate can only acquire the phase factor after the DW 
passage. 



^{xo^L)=e'^^e'^^{xo). 



(A3) 



The first factor represents the dynamical phase due to the 
spin precession in magnetic field, and the second factor is 
geometric Berry phase, independent of the DW structure 
in the adiabatic limit. Below we calculate the geometrical 
phase explicitly for the specific model of DW specified 
by Eqs. (|Aip and (|A2p . and compare it with the known 
theoretical results. 

We calculate transmission amplitudes t± for two 
spinors '^+(-) which are exact eigenstates of the Hamil- 
tonian for X < Xq. The spin in these states is polarized 
parallel ( ant i- parallel) to the magnetic field outside the 
DW, i.e. 



(V^±|a|^±)=± 



In the last equation we use the notation 



Vx 



2BxL/7TVf 



(A4) 



(A5) 



similar to Eq. (|6]). The transmission amplitudes are di- 
agonal elements of the inverse transfer matrix T defined 
in Eq. (dH), 



(A6) 



The unitary transfer matrix T is found following the same 
approach as in Sec. IV 



T = U exp {iTTCTx) exp {—iEL/vf) , 



(j = cos -V^2^fe-A^/^)2 

+^ / ^ , , An/ ^9 sm-Vr + fe-Ai9/7r)2. (A7) 



Substituting Eqs. (lATp and (IA4|) in Eq. (IA6|) we obtain 

.EL] 



t± =exp 



z Vf 



X exp 



^iiT 1 



Aer]x 



where we made an expansion 



(A8) 



(A9) 



valid in the adiabatic limit, r]^r]x ^ 1- Comparing 
Eq. (|A8l) with Eq. (|A3l) we identify the first factor as 
corresponding to the dynamical phase 



L/vp 



(AlO) 



due to the orbital motion and precession in the mag- 
netic field, and the second factor as the geometrical Berry 
phase. 



3 = ±7r(l — cos (Pb) , 



(All) 
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where (/^^ is an opening angle of the cone swept by the 
magnetic field in the DW. Equivalently, the phase S is 
half of the solid angle subtended by the closed contour 
in the B-space at the degeneracy point B = 0, in agree- 
ment with the general theory. We notice that the ge- 
ometrical phase is of opposite sign for two spin orien- 
tations tlj±. The resulting deviation of the precession 



angle from the one expected from Eq. (|AlQp has been 
observed experimentally for neutrons passing the region 
of the spiral magnetic field^^. We argue that similar phe- 
nomenon could be observed in graphene samples for car- 
riers passing the DW with the magnetic field rotation 
angle, = 27r. 
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